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Abstract. First results of simulations are presented which 
compute the dynamical evolution of a Chandrasekhar-mass 
white dwarf, consisting of equal amounts of carbon and oxy- 
gen, from the onset of violent thermonuclear burning, by means 
of a new two-dimensional numerical code. Since in the inte- 
rior of such a massive white dwarf nuclear burning progresses 
on microscopic scales as a sharp discontinuity, a so-called 
flamelet, which cannot be resolved by any numerical scheme, 
and since on macroscopic scales the burning front propagates 
due to turbulence, we make an attempt to model both effects 
explicitly in the framework of a finite-volume hydrodynam- 
ics code. Turbulence is included by a sub-grid model, follow- 
ing the spirit of large eddy simulations, and the well-localized 
burning front is treated by means of a level set, which allows 
us to compute the geometrical structure of the front more ac- 
curately than with previous methods. The only free parame- 
ters of our simulations are the location and the amount of nu- 
clear fuel that is ignited as an initial perturbation. We find that 
models in which explosive carbon burning is ignited at the cen- 
ter remain bound by the time the front reaches low densities, 
where we stopped the computations because our description of 
combustion is no longer applicable. In contrast, off-center ig- 
nition models give rise to explosions which, however, are still 
too weak for typical Type la supernovae. Possible reasons for 
this rather disappointing result are discussed. 

Key words: supernovae: general - physical data and processes: 
hydrodynamics - turbulence - nuclear reactions, nucleosynthe- 
sis, abundances - methods: numerical 



1. Introduction 

Despite the fact that considerable attention has recently been 
devoted to so-c alled sub-Chandrasekh armass models (see, e.g ., 
Ilben & Tutukov 199 it IWooslev & Weaver 1994 ILivne 19 96). 
the most popular progenitor model for at least some Type 
la supernovae is a massive white dwarf, consisting of car- 
bon and oxygen, which approaches the Chandrasekhar mass, 
Mch — 1.39 M Q , by a yet unknown mechanism, presumably 



accretion from a companion sta r, and is disrup ted by a ther- 
monuclear explosion (see, e.g., Woosley 1990 for a review). 
Arguments in favor of this hypothesis include the ability of 
these models to fit the observed spectra and light curves, in 
particular in versions that propagate the burning front mostly 
in form of a sub-sonic deflagration wave. 

However, not only is the evolution of massive white dwarfs 
to explosion very uncertain, leaving room for some diversity in 
the allowed set of initial conditions (such as the temperature 
profile at ignition), but also the physics of thermonuclear burn- 
ing in degenerate matter is complex and not well understood. 
The generally accepted scenario is that explosive carbon burn- 
ing is ignited either at the center of the star or off-center in a 
couple of ignition spots, depending on the details of the previ- 
ous evolution. For example, URCA-neutrino emission affects 
the temperature distribution as well as non-steady convection 
during the slow carbon burning phase prior to the explosion, 
and both effects are usually not treated in Type la models. Af- 
ter ignition, the flame is thought to propagate through the star 
as a sub-sonic deflagration wave which may or may not change 
into a detonation at low densities (around 10 7 g/cm 3 ). Numer- 
ical models with parameterized velocity of the burning front 
ha ve been very success ful, the prototype being the W7 model 
of INomoto et al. (19 84). However, these models do not solve 
the problem because attempts to determine the effective flame 
velocity from direct numerical simulations fail ed and gave 
velocities far too low for successful e xplosions (Livne 1993; 
IKhokhlov 1991 lArnett & Livne 1994 . This has led to some 
speculations about ways to change the de flagration into a su- 
personic detonation (Khokhl ov 1991albTcT) . 
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( Niemever & Hillebrandt 19951 
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Niemever & Wooslev 1997), we aim at a better understanding 
and a more accurate numerical treatment of thermonuclear 
deflagration waves in degenerate C+O white dwarfs, hoping 
to resolve the difficulties of this class of supernova models. In 
order to clarify the problem we shall first discuss in Sect. [2] 
the basic physics of deflagration fronts in the flamelet regime 
as well as attempts to model them by numerical simulations. 
In Sect. [3] we present a brief description of our new numerical 
method. The results of a series of computations are given 
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and discussed in Sect. |4] and, finally, Sect. [5] is devoted to a 
summary and outlook on future work. 



2. A model for turbulent combustion in the flamelet regime 

Nuclear burning in the degenerate matter of a dense C+O white 
dwarf, once ignited, is believed to propagate on microscopic 
scales as a conductive flame, wrinkled and stretched by local 
turbulence, but with essentially the laminar velocity. At high 
densities near the center of a Chandrasekhar-mass white dwarf 
the typical length scales for the width of the flame are a frac- 
tion of a millimeter. However, this is not the most relevant 
length scale. Due to the very high Reynolds numbers, which are 
of the order of 10 14 on typical macroscopic scales of 10 7 cm, 
macroscopic flows are highly turbulent and interact with the 
flame, in principle down to the Kolmogorov scale of 10~ 3 cm. 
This means that all kinds of hydrodynamic instabilities feed 
energy into a turbulent cascade, including the buoyancy- 
driven Rayleigh-Taylor instability and the s hear-driven Kelvin- 
Hel mholtz instability (see, e.g.. INiemever & Hillebrandt 19951 
and iHillebrandt & Niemever 19961) . Consequently, the picture 
that emerges is more that of a "flame brush" spread over the en- 
tire turbulent regime rather than a wrinkled flame surface. For 
such a flame brush, the relevant minimum length scale is the so- 
called Gibson scale, defined as the lower bound for the curva- 
ture radius of flame wrinkles caused by turbulent stress. Thus, 
if the thermal diffusion scale is much smaller than the Gibson 
scale (which is the case for the physical conditions of inter- 
est here) small segments of the flame surface are unaffected 
by large scale turbulence and behave as unperturbed laminar 
flames ("flamelets"). On the other hand side, since the Gibson 
scale is, at high densities, several orders of magnitude smaller 
than the integral scale set by the Rayleigh-Taylor eddies and 
many orders of magnitude larger than the thermal diffusion 
scale, both transport and burning times are determined by the 
eddy turnover times, and the effective velocity of the burning 
front is independent of the laminar burning velocity. This gen- 
eral picture has indee d been verified by laborato ry combustion 
experiments (see, e.g.. lAbdel-Gaved et al. 1987b . 

Nie mever (19951) and Nieme ver & Hillebrandt (19 951) have 
presented a numerical realization of this general concept. Their 
basic assumption was that wherever one finds turbulence this 
turbulence is fully developed and homogeneous, i.e. the turbu- 
lent velocity fluctuations on a length scale I are given by the 
Kolmogorov law 

/ I \ 1/3 

v(l)=v(L)\^-J , (1) 

where L is the integral scale, assumed to be equal to the 
Rayleigh-Taylor scale. Following the ideas outlined above, one 
can also assume that the thickness of the turbulent flame brush 
on the scale I is of the order of I itself. With these two as- 
sumptions and the definition of the Gibson scale one finds for 



^gibs r^l r^L — Art 

( I \ 1/3 

u(0^*(0-^(W> — (2) 

V 'gibs / 

and 

d t (l) ~ I, (3) 

where v(l g ib s ) — u\ am defines L^, u\. dm is the laminar burning 
speed and u t (l) is the turbulent flame velocity on the scale I. 

In a second step this model of turbulent combustion is cou- 
pled to a finite v olume hydro schem e such as the PPM-code 
PROMETHEUS jFrvxell et al. 1989b . Since in every finite vol- 
ume scheme scales smaller than the grid size cannot be re- 
solved, we express / g jb s in terms of the grid size A, the (un- 
resolved) turbulent kinetic energy q, and the laminar burning 
velocity: 

^= a U) • (4) 

Here q is determined from a sub-grid model Jdement 19931 
iNiemever & H illebrand t"l995l) and, finally, the effective turbu- 
lent velocity of the flame brush on scale A is given by 

u t {A) = max(wi am ,i;(A),i; flT ), (5) 

with v(A) = \f2q~ and vrt oc \JgA, where g is the local 
gravitational acceleration. 

The numerical scheme outlined so far has been applied to 
two-dimensional simulations of expl osive C-burning in Mch- 
mass C+O white dwarfs by INiemever & Hillebrandt (1 995) 
and Niemever et al. (1996) for a variety of initial conditions. 
They found that, with the exception of a model that was ignited 
off-center in a single small blob, all models were unbound at 
the end of the computations, and between 0.59 M Q (for central 
ignition) and 0.65 M of matter had been burnt, which would 
guarantee an observable, although weak, Type la supernova, 
but certainly not an average one, in particular because the ex- 
pansion velocities were too low. 

There are two possible interpretations of these results. Ei- 
ther they show that the "normal" Type la supernova is not the 
result of the explosion of a Mch-mass C+O white dwarf, caused 
by a deflagration wave, or the modeling of thermonuclear com- 
bustion as outlined so far is still insufficient. In this paper we 
investigate the second possibility and, in particular, present re- 
sults obtained by means of an improved numerical scheme. The 
main aim is to cure a problem that is common to all finite vol- 
ume schemes, namely that, although discontinuities are cap- 
tured automatically, they are smeared artificially over several 
grid zones. In many applications this broadening of disconti- 
nuities might not be a big problem but for modeling turbulent 
combustion one cannot be satisfied. 

The reasons for this are simple. First, as we discussed al- 
ready, it is important to reproduce the geometry of the flame 
front as accurately as possible, because its effective surface 
area also determines the rate of fuel consumption. If we think 
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of the flame in terms of a surface of a given fractal dimension 
the minimum scale determines essentially its effective velocity. 
Therefore it is important to resolve the small scales as accu- 
rately as possible. Secondly, since nuclear reactions are very 
sensitive to temperature and their rates depend on a very high 
power of T, mixing fuel and ashes numerically in certain mesh 
points on both sides of the front and smearing out the tempera- 
ture gradient there will lead to huge errors in the (local) energy 
generation rates. Therefore, if in a code one wants to follow 
the nuclear transmutations explicitly one is forced to guarantee 
that temperature jumps due to fast reactions are tracked with 
high precision. 

A possible way to reach both goals is to track the front by 
means of a level set function and to reconstruct the thermo- 
dynamical quantities in front and behind the the discontinu- 
ity from conservation laws. This approach has been utilized 
successfu lly for simulations of te rrestrial hydrogen combus- 
tion by ISmilianovski et al. (19971) and was reformulated and 
coded for the ap plication to thermonuc lear combustion in C+O 
white dwarfs by Reinec ke et al. (19981) . Here we apply a some- 
w hat simplified version of this c ode, i.e. we still assume, a s 
in lNiemever & Hillebrandt (19951) and lNiemever et al. (19 96). 
that the matter of the white dwarf is defined by two states, 
fuel and ashes, respectively, and that the transition from fuel 
to ashes is given by the instantaneous liberation of a certain 
amount of energy q once the material crosses the flame front. 
Thus the front-tracking by the level-set function can be inter- 
preted as a kind of passive advection without reconstruction. 
Nonetheless, in contrast to common front capturing schemes, 
our code is able to resolve complex structures of the flame on 
small scales. 

3. The numerical method 

The numerical code applied in this work is again based on 
the PPM code PROMETHEUS jFrvxell et al. 1989b combined 
with a sub-grid model for unresolved turbulence, as w as out- 
lined in Sect.|2l(see also Niemever & Hillebrandt (19951) 1. This 
includes that the propagation velocity of the burning front is 
obtained from the local turbulent kinetic energy as given by 
eq.0. 

As a model for the fla me, we use the so-called level set 
approach first described by Osh er & Sethian (19881) . The key 
concept of this model is to associate the flame front T(t) with 
the zero level set of a scalar function G(r,t): 



r(t) :={r\G(r,t) = 0} 



(6) 



This condition does not fully determine G. For numerical rea- 
sons it is appropriate to postulate also that G be positive in the 
burnt regions and negative before the front and that 



- The front is modeled as a discontinuity between fuel and 
ashes; this is an excellent approximation, if one compares 
the flame witdh and the typical cell size. A piecewise linear 
representation of the front is easy to obtain. 

- Complex topological situations, like the breaking of the 
burnt area into several disconnected regions, re-merging 
and highly convoluted front geometries are handled with- 
out problems. 

- The approach can be used for an arbitrary number of simu- 
lated dimensions. An extension to three-dimensional prob- 
lems is straightforward. 

G(r, 0) is defined by the initial conditions. During the sim- 
ulation, the front propagates with the velocity 



Aflame = V u + U t Tl 



(8) 



with respect to the underlying grid, where v u represents the ve- 
locity of the unburnt material and n is the front normal point- 
ing towards the unburnt material. This leads to the following 
expression for the temporal evolution of G in the vicinity of 
the front: 



dG 
~dt 



"Aflame VG = (v u Tl + U t )|VG| 



(9) 



To preserve the distance-function property, an additiona l re- 
initialization step is nece ssary; see lReinecke et al. (1998!) and 
Smiliano vski et al. (19971) for a detailed discussion. 

The level set propagation according to eq. [9] requires the 
knowledge of the velocity of the unburnt material v u at the 
front. Although this val ue can, in principle, be ca lculated in the 
cells cut by the front (Smi lianovski et al. 19971) . the extreme 
thermodynamical properties of the degenerate white dwarf ma- 
terial make this reconstruction very difficult. 

Fortunately, under the conditions typical for Type la super- 
novae, the velocity jump is quite small compared to the burning 
velocity, so that v u can be approximated by the average flow 
velocity with only small error. In this case, the hydrodynamical 
advection of G can be written in conservation form: 



L 



d( P G) 



v FP Gdf = 



(10) 



dv 



Therefore it is possible to treat the advection of the front caused 
by the fluid movement with the Godunov method built into 
PROMETHEUS. 

The propagation of the front, the species conversion and the 
energy release due to burning is done in an additional subrou- 
tine. In this step we apply 



IVGI = 1, 



(7) 



G' = G + Atu t JD 2 x 



D 2 



(11) 



which means that G is a signed distance function. 

This flame description is well suited for the combustion in 
Type la supernovae for a variety of reasons: 



where D x and D y represent discrete gradients. Next, the burnt 
volume fraction of all cells cut by the front are calculated, and 
the chemical composition and energy of these cells is adjusted. 
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Fig. 1. Initial front geometry for the centrally ignited models CI 
4. Results and discussion 

We have carried out numerical simulations in cylindrical rather 
than in spherical coordinates, mainly because it is much sim- 
pler to implement the level set on a Cartesian (r,z) grid. More- 
over, we save computer time since we avoid very small CFL 
time-steps which result in spherical coordinates near the center 
of the star if one wants to have good angular resolution. 

The grid we used maps the white dwarf onto 256x256 
mesh points, equally spaced for the innermost 226x226 
zones by A=1.5T0 6 cm, but increasing by 10% from zone 
to zone in the outer parts. The white dwarf, constructed in 
hydrostatic equilibrium for a realistic equation of state, has 
a central density of 2.9T0 9 g/cm 3 , a radius of 1.5T0 8 cm, 
and a mass of 2.8-10 33 g, id entical to the one used by 
iNiemever & Hillebrandt (19951) . The initial mass fractions of C 
and O are chosen to be equal, and the total binding energy turns 
out to be 5.4- 10 50 erg. Finally, in order to solve for the sub-grid 
turbulent kinetic energy, we have to specify its initial values on 
the grid and we assume a constant small value of q=10 12 erg/g, 
again in agreement with Niemever & Hillebrandt (1995). This 
latter assumption is justified because after a very short time q 
adjusts to the local turbulence, independent of its initial value. 
At low densities (p < 10 7 g/cm 3 ), the burning velocity of the 
front is set equal to zero because the flame enters the distributed 
regime and our physical model is no longer valid. However, 
since in reality significant amounts of mass may be burned at 
lower densities, the energy release obtained from our simula- 
tions is a lower limit only. On the other hand, it is difficult to 
estimate the error due to this simplification, mainly because the 
physics in the distributed regime of burning is not yet under- 
stood. It may happen that the flame just "stalls", in which case 
our numerical description would be appropriate. In the other 
extreme, even a deflagration-to-detonation transition (DDT) is 
not excluded which would change the energetics considerably. 

4.1. Centrally ignited models 

In this set of simulations we start the computations assuming 
that the inner 1.5-10 7 cm, corresponding to a mass of about 
O.O2M , are burnt. In addition, a sinusoidal perturbation is su- 



. t=0.00s, p c = 2.89e+09g/cm 3 
2.Q.10' r i 
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r[cm] 

andC3. 

perimposed producing one "finger" (model CI) or three fingers 
(model C3), respectively, as seeds for the Rayleigh-Taylor in- 
stability (see Fig.Q. 

Fig.|2]shows a sequence of snapshots of the front geometry 
and the velocity field for model CI (note that all scales change 
from snapshot to snapshot). The burnt matter is shaded and the 
front, given by the level set function, is plotted as a solid line. 
The evolution follows essentially one's intuition. First the "fin- 
ger" grows into a typical Rayleigh-Taylor structure, mainly in 
radial direction, giving rise to two vortices on both sides of 
the Rayleigh-Taylor head. There, the turbulence intensity in- 
creases and nuclear burning is accelerated, leading, in turn, to 
still faster motions. For a while, a burning blob even discon- 
nects from the much slower burnt matter near the center of the 
star. In part, this is due to our symmetry assumptions which 
force us to use reflecting boundary conditions on both coordi- 
nate axes, but it also shows that our code can handle multiply 
connected fronts. 

After about Is most of the matter of the white dwarf is 
expanding and roughly 0.3M Q have been burnt into Ni, too 
little to unbind the star. The structure of the burnt region be- 
comes very complex. Kelvin-Helmholtz whirls can be seen 
at the interfaces between rising hot matter and unburnt more 
slowly expanding fuel. Once the central density has dropped to 
310 7 g/cm 3 we stopped the calculations for the reason which 
was explained earlier. The amount of Ni produced then is 
nearly 0.33M Q . The obvious asymmetries in the flow patterns 
at late times are again a consequence of the imposed sym- 
metry conditions: Convective eddies are not bubbles but axi- 
symmetric rings, and while for an equidistant grid the cell vol- 
ume does not increase in z-direction, it does so in radial direc- 
tion proportional to r because of the axi-symmetry. 

Fig. |3] shows snapshots of model C3 where we started the 
computations with three "fingers" rather than one. The changes 
again are not unexpected. Due to the more complex initial per- 
turbation burning propagates faster in the beginning, mainly 
because the effective surface of now three Rayleigh-Taylor 
bubbles is bigger. The net effect, however, is that already af- 
ter about 0.6s the white dwarf expands and cools, and, in the 
end, the total amount of Ni produced is only marginally higher, 
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Fig. 2. Temporal evolution of front geometry and velocity distribution for model CI. The inequal spacing of the velocity arrows 
in the last snapshot results from the non-equidistant grid points. Note that all scales change from snapshot to snapshot. 
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Fig. 3. Temporal evolution of front geometry and velocity distribution for model C3. 
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Fig. 4. Initial front geometry for the off-center ignited models Bl and B5. 



by about 10%, as in model CI, but also this model remains 
weakly bound. 

4.2. Off-center ignition models 

In two more sim ulations we follow the work of 
iNiemever et al. (19961) and ignite the carbon and oxygen 
fuel off-center, first in one blob and then in five ones (models 
Bl and B5; see Fig. |4] for the initial configurations). Since 
these blobs can rise freely, not connected to the center by a 
symmetry condition, we expect that they move faster, thereby 
generating more turbulence and, consequently, also more rapid 
burning. Figs. [5] and |6] demonstrate that this expectation is 
indeed fulfilled by the simulations. In models Bl and B5, 
after 0.61s, peak velocities are typically higher than those of 
models CI and C3 by up to 50%. Because fuel is consumed 
very rapidly in Bl and B5 also the expansion sets in earlier. 
However, this does not extinguish the burning in the case of 
model B5. In contrast, because of the large surface of the 
various burning blobs fuel consumption progresses so fast 
that this model is unbound already after 0.5s, and continues 
to burn until the amount of Ni is 0.34M Q . Therefore, in our 
computations this model is the only one which gives rise to an 
observable supernova explosion. Model Bl, on the other hand 
side, ends also unbound, but its total energy at the end of the 
computations is close to zero. 

A final remark concerns numerical effects which are obvi- 
ously present in our simulations and which may affect the out- 
come. First, as was stated before, in three dimensions without 
artificial symmetry assumptions Rayleigh-Taylor eddies will 
move faster. Moreover, a multiple of rising blobs has a higher 
surface area than our Rayleigh-Taylor "rings". Both effects will 
increase the rate of fuel consumption and, therefore, will work 
in favor of more violent explosions. Moreover, due to finite 
numerical resolution, the minimum scale of eddies shown in 
the previous figures is certainly too large. For example, when- 
ever blobs disconnect from each other, this is because burn- 
ing strings in between them are not resolved. This, in turn 
means that our code still underestimates the total rate of nu- 
clear burning. Of course, we hope that to a certain extend the 
subgrid model for the unresolved turbulent scales takes care 





CI 


C3 


Bl 


B5 


Elapsed time [s] 


1.35 


1.43 


1.22 


1.16 


Total mass of 56 Ni [M©] 


0.23 


0.26 


0.27 


0.34 


Total energy [10 5() erg] 


-0.82 


-0.23 


0.04 


1.32 


Kinetic energy [10 erg] 


1.95 


2.00 


2.55 


3.11 



Table 1. Comparison of integral quantities at the end of the four 
simulations. The numbers given are only approximate because 
of some mass-loss from the computational grid at late stages. 



of these effects. However, since we do not resolve the small- 
est RT-unstable wavelengths the assumption of a Kolmogorov- 
spectrum on sub-grid scales may not be justified and, even 
more important, we may not have enough spatial resolution to 
reach the turbulent regime. In principle, this hypothesis could 
be tested by increasing the resolution significantly, making the 
simulations very expensive, and first results of simulations with 
A = 5 • 10 5 cm give (after t « 0.6s) indeed a considerably 
higher energy production rate than the models discussed above, 
indicating that the calculations are not fully converged. It ap- 
pears to be unlikely, however, that going to higher numerical 
resolution will change our main conclusions. 

A numerical comparison of all four models at the end of the 
simulation is given in Tabled 

4.3. A comparison with previous works 

The results outlined s o far should be compared to those 
of INiemever (1995b INiemever & Hillebrandt (19951) and 
INiemever et al. (19961) . because the main motivation behind 
our present work was to improve their numerical treatment of 
a turbulent burning front in the flamelet regime, hoping that 
such improvements would change their weak explosions into 
more violent ones. Now the contrary appears to be true. 

In order to make the comparison more quantitative we plot 
in Figs. through |9] the total energy, the nuclear energy gen- 
eration rate, and the turbulent sub-grid energy as functions of 
time for our models and for some models from the previous 
papers. Since with the exception of the front-tracking scheme 
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of the present computations both codes are identical, simula- 
tions with similar initial conditions allow us to compare the 
effects caused by the flame-propagatio n model. In this respect, 
our model Bl and the model IB of lNiemever (19951) . as well 
as B5 and his "slow 5B" are similar. The same holds for our 
model C3 and Niemeyer's CI, with the difference that the per- 
turbation wavelength in our simulation had to be larger due to 
limited angular resolution near the center of the star. 

The first major difference seen in Fig. is that, in con- 
trast to our expectations, the code with front-tracking in general 
produces less energy from nuclear burning than the other one. 
(Note that the energies plotted include already the energy liber- 
ated in the initial perturbation. Therefore the curves do not start 
from the same value at time t=Qs.) But a closer look at Figs. 
[8] and [9] re veals the reason. Fig. [8] gives the energy generation 
rate for the various models. In those models with front- tracking 
the energy generation rate rises more rapidly in the beginning, 
caused by the rapid increase of the turbulent sub-grid energy 
(see Fig. |9jl, in contrast to the models without front- tracking. 
This, in fact, is in agreement with ones expectations because 
the tracked front has more structure and, therefore, a larger sur- 
face area. However, already after 0.5s the total energy released 
by nuclear burning in all our present models is high enough 
to lead to an overall expansion of the white dwarf, as was dis- 
cussed in the previous subsections. Consequently, the sub-grid 
energy drops again and so does the nuclear energy generation 
rate. So in total less energy is produced in comparison to the 
models without tracking, which reach their peak values of both 
the energy generation and the sub-grid energy about a factor 
of two later and, therefore, expand later. An exception is our 
model B5. Here the energy generation rate rises so fast that the 
star is already unbound at the time when bulk expansion sets 
in. 

In addition to the higher sub-grid energy production during 
the early stages of the simulations, we observe that most of the 
turbulence is generated in a very thin region around the flame, 
resulting in extremely high turbulence intensities near the front 
and a further increase in the burning speed (Fig. HOt . This is 
caused by the fact that the transition between the rising bubbles 
of hot, light ashes and the dense, cold fuel is quite thin (ps 2 A), 
which leads to a well-localized shear flow and turbulent energy 
generation, as one would expect in reality. In Niemeyer's simu- 
lations, in contrast, the transition was smeared out over several 
cells because of the employed reactive-diffusive flame model; 
consequently, the turbulent flame propagation speed was un- 
derestimated. 

5. Summary and conclusions 

In the present paper we have presented first results obtained by 
applying a new numerical technique to thermonuclear explo- 
sions of Chandrasekhar-mass C+O white dwarfs. Our code dif- 
fers from those used previously in the atrophysics literature in 
that it follows the propagation of the burning front explicitly by 
means of a level set function. Our approach has the advantage 
that numerical diffusion of the front is largely avoided and that 




3.O10 7 6.O10 7 9.O10 7 1.2«10 ; 
r[cm] 



Fig. 10. Distribution of the turbulent velocity fluctuations 
(Model B5, after 0.4s) 

the structure of burning regions is resolved down to a few grid 
zones, which is an important issue in numerical simulations of 
combustion. Moreover, in principle, an extension of the code is 
possible, making use of the fact that the thermodynamic prop- 
erties in mixed cells can be reconstructed from the conservation 
laws of mass, energy, and momentum, if the burning velocity 
of the front is known. Such a code could also be applied to the 
Type la supernova problem and would, for the first time, allow 
to include individual nuclear reactions in a meaningful way, 
and we are presently working on this extension. 

Since our aim was to demonstrate that deflagrations in Mch 
white dwarfs can lead to explosions which have the properties 
of some (if not of typical) Type la supernovae, the outcome so 
far is disappointing. Our models produce even less nickel and 
less energy than those computed with simpler and less accurate 
numerical schemes. In retrospect the reason for this finding is 
obvious. For a "healthy" explosion it is not sufficient to accel- 
erate the burning front beyond what is predicted by previous 
numerical experiments, at least not in multi-dimensional mod- 
els. In fact, more rapid local burning may under certain circum- 
stances result in an expansion and cooling of the white dwarf 
before large amounts of nuclear fuel are consumed. Our cen- 
trally ignited models and the one with a single off-center blob 
are examples. 

The question, therefore, remains whether or no our simula- 
tions rule out the deflagration scenario for typical Type la su- 
pernovae. Possible short-comings of the present numerical ap- 
proach, such as the artificial symmetry assumption or the still 
insufficient numerical resolution, have already been mentioned. 
Here improvements will be possible with increasing computer 
power in the near future. It might also be worthwhile to inves- 
tigate the behaviour of the simulations with varied parameter 
sets, e.g. for different chemical compositions or a slightly in- 
creased turbulent flame speed, since both of these quantities 
are not known exactly. 

A second open question concerns our model of turbulent 
combustion. Here it is not clear at all, if one of our basic as- 
sumptions, namely that in the presence of reactions the turbu- 
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lence spectrum is given by the Kolmogorov law remains valid. 
Although in the limit of very high turbulence intensity experi- 
ments seem to support our hypothesis more work needs to be 
done. Also, significant burning is still possible at low densities 
in the so-called distributed regi me, and even a transition to a 
detonation is not ruled out there ( Niemever & Wooslev 1997h . 

Finally, white dwarfs at the onset of the explosion might 
look rather different from the ones we have used as initial con- 
ditions. URCA-neutrino emission and non stationary convec- 
tion during the evolution just prior to the explosion have al- 
ready been mentioned as likely sources of large uncertainties. 
Also, whether the star reaches the Chandrasekhar-mass by ac- 
cretion or by merging with a companion will make a big differ- 
ence. For example, rotation in some form may directly affect 
the propagation of the deflagration front. All these questions 
need to be investigated before we can safely dismiss the defla- 
gration models. 
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